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■ Abstract. Effective Field Theory (EFT) provides a powerful framework that exploits 
a separation of scales in physical systems to perform systematically improvable, 
modcl-indcpcndcnt calculations. Particularly interesting are few-body systems with 

■ short-range interactions and large two-body scattering length. Such systems display 
£NJ \ remarkable universal features. In systems with more than two particles, a three-body 

force with limit cycle behavior is required for consistent renormalization already at 
leading order. We will review this EFT and some of its applications in the physics 
of cold atoms and nuclear physics. In particular, we will discuss the possibility of an 
^sO • infrared limit cycle in QCD. Recent extensions of the EFT approach to the four-body 

system and iV-boson droplets in two spatial dimensions will also be addressed. 

o 
in 
o _ 

1. Introduction 

The Effective Field Theory (EFT) approach provides a powerful framework that exploits 
the separation of scales in physical systems [T|. Only low-energy (or long-range) degrees 
of freedom are included explicitly, while all others are parametrized in terms of the 
most general local contact interactions. This procedure exploits the fact that a low- 
energy probe of momentum k cannot resolve structures on scales smaller than 1/k. 
(Note that H = 1 in this talk.) Using renormalization, the influence of short-distance 
physics on low-energy observables is captured in a few low-energy constants. Thus, the 
EFT describes universal low-energy physics independent of detailed assumptions about 
the short-distance dynamics. All physical observables can be described in a controlled 
expansion in powers of kl, where / is the characteristic low-energy length scale of the 
system. The size of / depends on the system under consideration: for a finite range 
potential, e.g., it is given by the range of the potential. For the systems discussed here, 
I is of the order of the effective range r e . 

We will focus on applications of EFT to few-body systems with large S-wave 
scattering length \a\ ^> I. For a generic system, the scattering length is of the same 
order of magnitude as the low-energy length scale I. Only a very specific choice of 
the parameters in the underlying theory (a so-called fine tuning) will generate a large 
scattering length a. Nevertheless, systems with large a can be found in many areas 
of physics. The fine tuning can be accidental or it can be controlled by an external 
parameter. Examples with an accidental fine tuning are the S-wave scattering of 
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nucleons or of 4 He atoms. The scattering length of alkali atoms close to a Feshbach 
resonance can be tuned experimentally by adjusting the external magnetic field. At very 
low energies these systems behave similarly and show universal properties associated 
with the large a [2 . We will start with a brief review of the EFT for few-body systems 
with large a and then discuss some applications in nuclear and atomic physics. 

2. Three-body system with large scattering length 

For typical momenta k ~ 1/1 a l> the EFT expansion is in powers of l/\a\; in leading 
order we can set / = 0. We start with a two-body system of non-relativistic bosons 
with large S-wave scattering length a and mass m. The generalization to fermions is 
straightforward. In the following, we will refer to the bosons simply as atoms. At 
sufficiently low energies, the most general Lagrangian may be written as: 



where the dots represent higher-order derivative terms which are suppressed at low- 
energies. The strength of the two-body interaction Cq is determined by the scattering 
length a, while Dq depends on a three-body parameter to be introduced below. 
For momenta k of the order of the inverse scattering length l/|a|, the problem is 
nonperturbative in ka. The exact two-body scattering amplitude can be obtained 
analytically by summing the so-called bubble diagrams with the Cq interaction term. 
The Dq term does not contribute to two-body observables. After renormalization, 
the resulting amplitude reproduces the leading order of the well-known effective range 
expansion for the atom-atom scattering amplitude: fAA(k) = (— 1/a — ik)^ 1 , where the 
total energy is E = k 2 /m. If a > 0, Jaa has a pole at k = i/a corresponding to a shallow 
dimer with binding energy B 2 = l/(ma 2 ). Higher-order derivative interactions are 
perturbative and give the momentum- dependent terms in the effective range expansion. 

We now turn to the three-body system. Here, it is useful to introduce an auxiliary 
field for the two-atom state (see Ref. j2] for details). At leading order, the atom- 
dimer scattering amplitude is given by the integral equation shown in Fig. A solid 



line indicates a single atom and a double line indicates the interacting two-atom state 
(including rescattering corrections). The integral equation contains contributions from 
both the two-body and the three-body interaction terms. The inhomogeneous term is 
given by the first two diagrams on the right-hand side: the one-atom exchange diagram 




(1) 




Figure 1. Integral equation for the atom-dimer scattering amplitude. Single (double) 
line indicates single atom (two-atom state). 
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and the three-body interaction. The integral equation simply sums these diagrams to 
all orders. After projecting onto S-waves, we obtain the equation 

3a 71 Jo _i/ a+ ^3q 2 /4-mE - it 

for the off-shell atom-dimer scattering amplitude with the inhomogeneous term 

1 , / p 2 +pk + k 2 -mE\ H(A) 

The logarithmic term is the S-wave projected one-atom exchange, while the term 
proportional to H(A) comes from the three-body force. The S-wave atom-dimer 
scattering amplitude /ad(&) = [kcotS — ik}~ 1 is given by the solution T evaluated 
at the on-shell point: fAD{k) = T(k, k; E) where mE = 3k 2 /4 — 1/a 2 . The three-body 
binding energies B 3 are given by those values of E < for which the homogeneous 
version of Eq. (J2J) has a nontrivial solution. 

Note that an ultraviolet cutoff A has been introduced in (J2J). This cutoff is required 
to insure that Eq. (J2J) has a unique solution. All physical observables, however, must be 
independent of A, which determines the behavior of if as a function of A [3]: 

_ cos[goln(A/A„) + arctansp] 
cos[sohi(A/A*) — arctanso] ' 
where So = 1.00624 is a transcendental number and A* is a three-body parameter 
introduced by dimensional transmutation. This parameter cannot be predicted by the 
EFT and must be determined from a three-body observable. Note also that H(A) is 
periodic and runs on a limit cycle. When A is increased by a factor of exp(7r/ s ) ~ 22.7, 
H(A) returns to its original value. 

In summary, two parameters are required to specify a three-body system at leading 
order in l/\a\: they may be chosen as the scattering length a (or equivalently B 2 if 
a > 0) and the three-body parameter A* 0. 



3. Universal properties of few-body systems with large scattering length 

This EFT confirms and extends the universal predictions for the three-body system first 
derived by Efimov [4 . The best known example is the Efimov effect, the accumulation 
of infinitely many three-body bound states at threshold as a — > ±oo. Universality also 
constrains three-body scattering observables. The atom-dimer scattering length, e.g., 
can be expressed in terms of a and A* as jH Ej 

ai2 = a (1.46 -2.15 tan[s ln(aA») +0.09]) (1 + 0(1 / a)) , a>0. (5) 

Note that the log-periodic dependence of H on A* is not an artifact of the renorma- 
lization and also shows up in observables like au- This dependence could, e.g., be 
tested experimentally with atoms close to a Feshbach resonance by varying a. Similar 
expressions can be obtained for other three-body observables such as scattering phase 
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shifts as well as rates for three-body recombination and dimer relaxation (including the 
effects of deeply-bound two-body states) 0. 

Since up to corrections of order l/\a\, low-energy three-body observables depend 
on a and A* only, they obey non-trivial scaling relations. If dimensionless combinations 
of such observables are plotted against each other, they must fall close to a line 
parametrized by A* j^JU]. An example of a scaling function relating 4 He trimer ground 
and excited state energies B% and S3 is shown in the left panel of Fig. [21 (A related 




B 3 /B 2 a 12 /a. 



Figure 2. Scaling functions relating the 4 He trimer ground and excited state energies 
(left panel) and the 4 He trimer ground state and the atom-dimer scattering length 
(right panel). The data points are calculations using various methods and 4 He 
potentials (see Ref. and references therein). Note that 0,3 = l/y/mBz. 

scaling function was obtained in Ref. [H].) The data points show calculations using 
various approaches and 4 He potentials. Since these potentials have approximately the 
same scattering length but include different short-distance physics, the points on this 
line correspond to different values of A*. The small deviations of the potential model 
calculations are mainly due to effective range effects. They are of the order r e /a rs 10% 
and can be calculated at next-to-leading order in EFT. The calculation corresponding to 
the data point far off the universal curve can easily be identified as problematic since the 
deviation from universality by far exceeds the expected 10%. The right panel shows the 
scaling function relating the 4 He trimer ground state energy and the atom-dimer 
scattering length au- A similar scaling relation is observed in nuclear physics between 
the spin-doublet neutron-deuteron scattering length and the triton binding energy and 
is known as the Phillips line. 

Recently, we have extended the effective theory for large scattering length to the 
four-body system [HJ. It is advantageous in this case to use an effective quantum 
mechanics framework, since one can start directly from the well-known Yakubovsky 
equations. We have shown that no four-body parameter enters at leading order in 
l/\a\. Therefore renormalization of the three-body system automatically guarantees 
renormalization of the four-body system. Consequently, there are universal scaling 
functions relating three- and four-body observables as well. As an example, we display 
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the correlation between the trimer and tetramer binding energies in Fig. El The left panel 
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Figure 3. Scaling function relating the 4 He trimer and tetramer ground state (left 
panel) and excited state (right panel) energies. The data points are calculations using 
various methods and 4 He potentials. 

shows the correlation between the ground state energies while the right panel shows the 
correlation between the excited state energies. The data points are calculations using 
various methods and 4 He potentials (see Ref. [H] and references therein). We conclude 
that universality is well satisfied in the three- and four-body systems of 4 He atoms. 

The existence of these scaling functions is a universal feature of systems with 
large scattering length and is independent of the details of the short-distance physics. 
Similar correlations between few-body observables appear for example also in nuclear 
systems. The correlation between the triton and a particle binding energies (the Tjon 
line) can be explained using a similar effective theory ^H]- The spin-singlet and spin- 
triplet scattering lengths a s and a t for nucleons are both significantly larger than the 
range of the nuclear force. This observation can be used as the basis for an EFT 
approach to the few-nucleon problem in which the nuclear forces are approximated 
by contact interactions with strengths adjusted to reproduce the scattering lengths a s 
and a t [TH [TJ|. This EFT does not contain explicit pion degrees of freedom and is 
sometimes referred to as pionless EFT. The EFT involves an isospin doublet N of 
Pauli fields with two independent 2-body contact interactions: N^o~iN c N c ^aiN and 
NWkN c N°WkN, where iV c = a 2 r 2 N*. Renormalization in the 2-body sector requires 
the two coupling constants to be adjusted as functions of A to obtain the correct values 
of a s and a t . Renormalization in the 3-body sector requires a 3-body contact interaction 
AH o~iN c N°i 'ajNN' o~iO~jN with a coupling constant proportional to (J1J) [12J. Thus the 
renormalization involves an ultraviolet limit cycle. The scaling-violation parameter A* 
can be determined by using the triton binding energy B t as input. Effective range effects 
and other higher order corrections can be treated in perturbation theory [TH] . 

In Fig. we show the result for the Tjon line with a s and as input (lower line) 
and a s and at as two-body input (upper line). Both lines generate a band that gives a 
naive estimate of higher order corrections in £/\a\. We also show some calculations using 
phenomenological potentials fl] and a chiral EFT potential with explicit pions [15 j . All 
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Figure 4. The correlation between the binding energies of the triton and the a- 
particle (the Tjon line). The lower (upper) line shows our leading order result using 
a s and Bd (a s and at) as two-body input. The gray circles and triangles show 
various calculations using phcnomcnological potentials without or including three- 
nucleon forces, respectively. The squares show the results of chiral EFT at NLO 
for different cutoffs while the diamond shows the N 2 LO result. The cross shows the 
experimental point. 

calculations with interactions that give a large scattering length should lie close to this 
line. Different short- distance physics and/or cutoff dependence should only move the 
results along the Tjon line. This can for example be observed in the NLO results with 
the chiral potential indicated by the squares in Fig. HJ 

4. An infrared renormalization group limit cycle in QCD 

The low-energy few-nucleon problem can also be described by an EFT that includes 
explicit pion fields. Such an EFT has been used to extrapolate nuclear forces to the 
chiral limit of QCD in which the pion is exactly massless ^Hl E|- The extrapolation 
to larger values of m n predicts that at diverges and the deuteron becomes unbound at 
a critical value in the range 170 MeV < m n < 210 MeV. It is also predicted that a s is 
likely to diverge and the spin-singlet deuteron become bound at some critical value of 
not much larger than 150 MeV. Both critical values are close to the physical value 
m n = 135 MeV. 

Chiral extrapolations can also be calculated using the EFT without pions JB]- The 
inputs required are the chiral extrapolations a s (m,r), a t (m n ), and B t (m n ), which can 
be calculated using an EFT with pions. As an illustration, we take the central values 
of the error bands for the inverse scattering lengths l/a s (m n ) and l/a t (m n ) from the 
chiral extrapolation in Ref. jEj. Since the chiral extrapolation of the triton binding 
energy B t (m n ) has not yet been calculated and since A* should vary smoothly with m n , 
we approximate it by its physical value A* = 189 MeV for = 138 MeV. In Fig. EJ 
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Figure 5. The binding momenta n — (771B3) 1 / 2 of pnn bound states as a function 
of TO^f. The circles are the triton ground state and excited state. The crosses are 
the physical binding energies of the deuteron and triton. The dashed lines are the 
thresholds for decay into a nucleon plus a deuteron (left curve) or a spin-singlet 
deuteron (right curve). 



we show the resulting 3-body spectrum in the triton channel as a function of m^. Near 
« 175 MeV where the decay threshold comes closest to k = 0, an excited state of the 
triton appears. This excited state is a hint that the system is very close to an infrared 
limit cycle. In the case illustrated by Fig. the value of m T at which a t diverges is 
larger than that at which a s diverges. If they both diverged at the same value of m n , 
there would be an exact infrared limit cycle. 

We conjecture that QCD can be tuned to this infrared limit cycle by adjusting the 
up and down quark masses m u and m d [THj. As illustrated in Fig. the tuning of m n , 
which corresponds to m u + m<2, is likely to bring the system close enough to the infrared 
limit cycle for the triton to have one excited state. We conjecture that by adjusting 
the two parameters m u and to critical values, one can make a t and a s diverge 
simultaneously. At this critical point, the deuteron and spin-singlet deuteron would both 
have zero binding energy and the triton would have infinitely-many increasingly-shallow 
excited states. The ratio of the binding energies of successively shallower states would 
rapidly approach a constant Aq close to 515. It may be possible to use a combination 
of lattice gauge theory and EFT to demonstrate the existence of this infrared RG limit 
cycle in QCD. 

5. iV-boson droplets in 2D 

In this section, we consider the universal properties of weakly interacting bosons with 
large scattering length (or equivalently a shallow dimer state) in two spatial dimensions 
(2D) [TH]. In particular, we consider self-bound droplets of 1) bosons interacting 

weakly via an attractive, short-ranged pair potential. Our analysis relies strongly on the 
property of asymptotic freedom of 2D bosons with an attractive interaction. 
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In 2D, any attractive potential has at least one bound state. For the potential 
— gS 2 (r) with small g, there is exactly one bound state with an exponentially small 
binding energy, B 2 ~ A 2 exp (— Air/g) , where A is the ultraviolet momentum cutoff 
(which is the inverse of the range of the potential). Asymptotic freedom provides an 
elegant way to understand this result. In 2D non-relativistic theory, the four-boson 
interaction g(ijj^ijj) 2 is marginal. The coupling runs logarithmically with the length 
scale R, and the running can be found by performing the standard renormalization 
group (RG) procedure. For g > 0, the coupling grows in the infrared, in a manner 
similar to the QCD coupling. The dependence of the coupling on the length scale R is 
given by 

(6) 

so the coupling becomes large when R is comparable to the size of the two body 
bound state B^ 1 ^ 2 . This is in essence the phenomenon of dimensional transmutation: a 
dynamical scale is generated by the coupling constant and the cutoff scale. It is natural, 
then, that B2 is the only physical energy scale in the problem: the binding energy of 
three-particle, four-particle,... bound states are proportional to B 2 . In contrast to three 
spatial dimensions, there is no Efimov effect (or Thomas collapse) and no three-body 
parameter is required. The A^-particle binding energy B^, however, can be very different 
from B 2 if A" is parametrically large. We use the variational method to calculate the size 
of the bound state. For a cluster of a large number of bosons, one can apply classical 
field theory. We thus have to minimize the expectation value of the Hamiltonian with 
respect to all field configurations ip{r) satisfying the constraint A" = fd 2 rip^ifj . This is 
equivalent to a Hartree calculation with the running coupling constant g(R) instead of 
the bare one. In the limit of a large number A" of particles in the droplet, some exact 
predictions can be obtained [TH] . 

The system possesses surprising universal properties. Namely, if one denotes the 
size of the A^-body droplet as Rn, then at large A" and in the limit of zero range of the 
interaction potential [T9] : 

fljv+i/fljv « 0.3417, B N+1 /B N « 8.567, N > 1 . (7) 

The size of the bound state decreases exponentially with A^: adding a boson into an 
existing A^-boson droplet reduces the size of the droplet by almost a factor of three. 
Correspondingly, the binding energy of A" bosons B^ increases exponentially with N. 
This implies that the energy required to remove one particle from a A^-body bound 
state (the analog of the nucleon separation energy for nuclei) is about 88% of the total 
binding energy. This is in contrast to most other physical systems, where separating 
one particle costs much less energy than the total binding energy, provided the number 
of particles in the bound state is large. The 1/A^-corrections to Eqs. (J7J) are calculable. 

For the universal predictions (|IJ) to apply in realistic systems with finite-range 
interactions, the A^-body bound states need to be sufficiently shallow and hence have 
a size Rn large compared to the natural low-energy length scale I. Depending on the 
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physical system, I can be the van der Waals length l v dw-, the range of the potential, or 
some other scale. As a consequence, Eqs. (J7J) are valid in such systems for N large, but 
below a critical value, 

1 « N < iV cri t ~ 0.931 hx{R 2 /l) + 0(1) . (8) 

At N = N CTit the size of the droplet is comparable to I and universality is lost. If there 
is an exponentially large separation between R 2 and /, then A cr i t is much larger than 
one and the condition (JHJ) can be satisfied. 

We can compare our prediction with exact few-body calculations for N = 3,4. 
While the 1/A-corrections are expected to be relatively large in this case, we can 
estimate how the universal result for Bjq + i/B^ is approached. The three-body system 
for a zero-range potential in 2D has exactly two bound states: the ground state with 
^ 0) = 16.522688(1)52 and one excited state with = 1.2704091(1)52 |2m ED E2| ■ 
Similarly, the four-body system for a zero-range potential in 2D has two bound states: 
the ground state with B^ = 197.3(1)S 2 and one excited state with B^ = 25.5(1)_B 2 
[22]. The prediction applies to the ground state energies and Bf\ The ratio 
-B3 / B 2 ~ 16.5 is almost twice as large as the asymptotic value (JZj), while the ratio 
B^/B^ pa 11.9 is already considerably closer. These deviations are expected for such 
small values of N. Note, however, that the ratio of the root mean square radii of the 
two- and three-body wave functions is 0.306 [2T], close to the asymptotic value (J7J). 

It would be interesting to test the universal predictions ((Z|) both theoretically and 
experimentally for N > 4. On the theoretical side, Monte Carlo techniques appear to 
be a promising avenue. Furthermore, the experimental realizability of self-bound 2D 
boson systems with weak interactions should be investigated. According to the analysis 
of Ref. [221, the 1/N corrections to Eqs. (J7J) are small for iV > 6. Using (jHJ), this requires 
Rili 3> 600. We are not aware of any physical system that satisfies this constraint. 
However, such a system could possibly be realized close to a Feshbach resonance where 
i?2 can be made arbitrarily large. 

6. Conclusions 

We have discussed the EFT for few-body systems with short-range interactions and 
large scattering length a. The renormalization of the three-body system with large a 
in three spatial dimensions requires a one-parameter three-body force governed by a 
limit cycle already at leading order in the expansion in l/\a\. As a consequence, two 
parameters are required to specify a three-body system: the scattering length a (or the 
dimer binding energy B 2 ) and the three-body parameter A*. Once these two parameters 
are given, the properties of the three- and four-body systems are fully determined at 
leading order in l/\a\. 

The large scattering length leads to universal properties independent of the short- 
distance dynamics. In particular, we have discussed universal expressions for three-body 
observables and universal scaling functions relating various few-body observables for 
both atomic and nuclear systems. A more detailed account can be found in Ref. |2j. 
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The success of the pionless EFT demonstrates that physical QCD is close to an 
infrared limit cycle. We have conjectured, that QCD could be tuned to the critical 
trajectory for the limit cycle by adjusting the up and down quark masses. The limit 
cycle would them be manifest in the Efimov effect for the triton [THj . 

In two spatial dimensions, the three-body parameter A* does not enter at leading 
order in the expansion in l/\a\ and iV-body binding energies only depend on B 2 . The 
asymptotic freedom of non-relativistic bosons with attractive interactions in 2D leads 
to some remarkable universal properties of iV-body droplets 

The three-body effects discussed here will also become relevant in Fermi systems 
with three or more spin states. Future challenges include universality in the iV-body 
problem for N > 4, effective range corrections, and a large scattering length in higher 
partial waves. (See Ref. [2j for more details.) A large P-wave scattering length appears 
for example in nuclear Halo systems such as 6 He |23| . 
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